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Abstract. Explicit formulae for the Bezier coefficients of the constrained dual Bern- 
stein basis polynomials are derived in terms of the Hahn orthogonal polynomials. 
Using difference properties of the latter polynomials, efficient recursive scheme is ob- 
tained to compute these coefficients. Applications of this result to some problems of 
CAGD is discussed. 
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1. Introduction 

(k I) 

Let IL\ ' , where k + l < n, be the space of all polynomials of degree < n, whose derivatives 



of order < k — 1 at t = 0, as well as derivatives of order < I — 1 at t = 1, vanish. In particular, 
II n := IL?' ^ is the space of the unconstrained polynomials of degree at most n. Obviously, 
dim = n — k — I + 1, and the Bernstein polynomials {-£>£■, B^ +l , . . . , -B"_ ; } , where 



B?(x) = f - z) n-i (0 < % < n), 

form a basis of this space. There is a unique dual constrained Bernstein basis of degree n, 
{D^ l \x-a,^ i)W\x;a,a..,I)<"f( I ;«,ffl}cnW), 

satisfying the relation /x)| n ' fc '^ ) = (i,j = k,k + l,...,n — l), where 5ij equals 1 if 
i = j and otherwise, and the inner product (•, •} is given by 

(1.1) (f,g):= I' (l-x) a xPf(x)g(x)dx (a, /3 > -1). 



o 



The unconstrained dual basis Df(x; a, f3) corresponds to the case k = I = 0, hence D™(x; a, /3) = 

Dual Bernstein polynomials D ( f l ' k ' l \x; a, (3) were studied by Ciesielski [1] (case a = (3 = 
and k = I = 0), also by Jiittler [6] (case a = f3 = and general k = I), Rababah and 
Al-Natour [11] (k = I = and general a, (3) and [12] (general a, /3 and k = I). 

In [8], we have given explicit expressions for D™(x; a, (3) in terms of the shifted Jacobi 

polynomials R^ a '^\x), which, as it is well known, form a set of orthogonal polynomials with 
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respect to the inner product (jl.ip . See also (|2,3p . In [IB], we show that the constrained 

polynomials D^' (x; a, (3) are related to the unconstrained dual Bernstein polynomials of 
degree n—k — l, with parameters a + 21 and (3 + 2k. Consequently, using the above-mentioned 
results of [8], we obtain an explicit orthogonal representation, as well as degree elevation 
formula for the constrained dual Bernstein basis polynomials. The latter result plays a crucial 
role in the algorithm of multi-degree reduction of Bezier curves, given in the cited paper. 

In some computer-aided geometric design applications, like curve intersection using Bezier 
clipping (see, e.g., [HJ [2j [10]), polynomial approximation of rational Bezier curves [9], or 
degree reduction of Bezier curves (see, e.g., [16J and the papers cited therein), it is desirable 
to have the Bezier-Bernstein representation of the dual constrained polynomials, 

n—l 

D^' k ' l \x; a,P)=J2 Cijin, k,l,a, (3) B?(x). 

j=k 

In this paper, we give an efficient algorithm to compute the coefficients Cij(n,k,l,a, f3) for 
i,j = 1,2... , n, which has 0(n 2 ) complexity. This should be compared to a rather compli- 
cated recursive algorithm given in [BJ and its generalization of |12j . 
In the sequel, we use the notation 

fe-i 

(c)k:=]J(c + j) (k>0). 
3=0 

The generalized hypergeometric series is defined by (see, e.g., [TJ §2.1], or [7, §1-4]) 



r ± s 



oo 

ai, . . . , a r 



= (ai)fc ■ ■ ■ (a r )k k 
■'■ ^ kl(bx)k . . . (b s )k ' 



h,...,b 

where r, s G Z + and oi, . . . , a r , b\, . . . , b s , z G C. 

2. Dual Bernstein polynomials 
Bernstein polynomial basis of degree n (n G N) is given by 

B»(ar) : = ( n \ x i (1 - x) n ^ (0 < i < n). 



Recall that 



(2.1) R^(x) := {a+ J )m 2 F X ( " m ' m + + a + P + 1 



1 — x 



are the shifted Jacobi polynomials \13\ p. 280, Eq. (46)]. For a > — 1 and f3 > —1, polynomials 
(|2.ip are orthogonal with respect to the scalar product (jl.ip (see, e.g., [131 P- 273]); more 
specifically, 

SW>) =6 mn B(a + l,P + 1) ^%^7% K n > 0) 
/ m\{2m/a + l){a) m 

with a := a + (3 + 1 and i?(A, /u) = r(A)T(/i)/r(A + (i). Recall that the Hahn polynomials are 
defined by §9.5] 

,„m n / a An f -m,m + a + P + l,-x 

(2.2) Q m (x;a,P,N) := 3 F 2 I ^ + 1 _^ 



1 (m = 0, 1,... ,JV; N G N). 
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Bernstein polynomials have the following representation in terms of shifted Jacobi polynomial 
basis: 

(2-3) Bf{x) = ( n ) (a + 1)^08 + 1), £ J+nffffi & n ^ 4* 

where < i < n. This result follows from a formula given in p3], using [U Cor. 3.3.5]; see 
also Eq. (5.4)]. Conversely, shifted Jacobi polynomials have the following representation 
in terms of Bernstein basis [3]: 

(2.4) R^\ x )=^±^i f2Qi(n-j;aJ,n) B]{x) (0 < i < n). 
Associated with the basis Bf(x) there is a unique dual Bernstein polynomial basis of degree 

n, 

D%(x; a, P), D?(x; a, /3), . . . , D%(x; a, /?), 

defined so that 

(2.5) (D?,B?) = 6ij (i,j = 0,l,...,n). 

In [8], we gave a number of properties of polynomials Df(x; a, /3), including the expansion 

- Bia 1) B-iy (2J {: : q^. «■ »)*r>w. 

where < i < n, as well as the following alternative formula representing the dual Bernstein 
polynomial Df(x;a, (3) (0 < i < n) as a "short" linear combination of Jacobi polynomials 
with shifted parameters: 

(27) D n (x-aB)- (-!)""> + l)n A j-i), R («,f>+j+i) , , 

In CAGD applications, it is desirable to express the dual basis polynomials in terms of the 
Bernstein basis. We have the following result. 

Theorem 2.1. For i = 0, 1, . . . , n, we have 

n 

(2.8) D?(x;a,P) = £ cyfa, a, 0) B?(x), 

j=0 

where the Bezier coefficients Cij(n,a,f3) (0 < j < n) are given by any of the following two 
forms: 

(2.9) 

m=0 

(2.10) 

/ (-l) n (a + l) n (-a-n)i >A (-i) ft . 

C *>' a ' # := £(a + l,/3 + l)n!(/3 + l) t g (-a - n) h Qn ~ h{n ~ r,*,P + h + l,n). 
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Proof. Formula f|2.8j) with Cij(n,a,f3) given by (|2,9p (respectively, (|2.10p ) follows by inserting 
expansion (12. 4j) in f)2.6f) (respectively, (|2.7|l ) and doing some algebra. □ 

Now, we show that the coefficients Cij(n,a, f3) = 1,2... ,n) given in Theorem 12.11 can 
be computed recursively, with a cost 0(n 2 ). 

Theorem 2.2. Quantities Cij = Cij(n,a,f3) defined in (12. 9|) . (I2.10P satisfy the following 
recurrence relation: 



-'J 



(2.H) + aj-i + A(j) c iJ+1 - B(i) Ci_ij}, 

where < i < n — 1, < j < n and 

A(/i) := (/i - n)(/i + /3 + 1), := h(h - n - a - 1). 

T/ie starting values are 

(0 , s (a + i) n (/3 + 2) n 

(2 - 12) ^ = n!S(a + l,/J + l) ' (a + !)„_,■(/? + 2), ( ° ^ n >" 



Proof. Observe that the expression on the right-hand side of (12. 9p is the linear combination 
of products Q m (i; (3, a, n) Q m (j; /3, a, n), with the coefficients independent of i and j. Hence, 
the proof of the Theorem can be based on an idea similar to the one used in |16l Thm A6] to 
prove a recursion satisfied by the dual discrete Bernstein polynomials. 

First, recall that the Hahn polynomials satisfy the difference equation (see, e.g., §9.5]), 

£x Q m (x; f3, a, N) = m(m + a) Q m (x; (3, a, N), 

where the operator £^ is given by 

£x v( x ) = un{x) y(x + 1) - c N (x) y(x) + b N (x) y(x - 1) 

with 

a N (x) := (x - N)(x + /3 + 1), b N (x) := x(x - a - N - 1), 
Cjv(x) := a N (x) + b N (x). 
Now, using equation (|2.9|) we easily obtain that 

(2m/<r + l)(/3 + l) 



n B(a + l, p + l)m\(a + l) 

m=0 ' 



x m(m + a)Q m (i;(3,a,n) Q m (j; p,a,ri) 



C n r- ■ 



This implies equation (|2.1ip . 

By (I2.10P and (12. 2p . the sum in coj reduces to z^l(J — n, n + cr + 1; a + 1; 1) which can be 
summed by the Chu-Vandermonde formula (see, e.g., [7, (1.5.4)]). Hence follows (12.12p . □ 
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Remark 2.3. Rababah and Al-Natour |12j have given the formula 
Ci ,(n,a,/3)=(-l)^'Q ^ 1 

min(i ,j) . , \ / r \ — 1 

EJ n + a + n \ n + a — h\ 
'n + /3 + /i + l\/n + a — /i 



n — m J \n + a — m 
which contains the earlier result of Jiittler [6] for a = /3 = 0. 

3. Constrained dual Bernstein polynomials 

Let Iln , where /c + / < n, be the space of all polynomials of degree < n, whose derivatives 
of order < k — 1 at t = 0, as well as derivatives of order < / — 1 at t = 1, vanish: 

IlW) : = |P € n n : PW(0) = (0 < * < Jfe - 1); P (j) (l) =0 (0 < j < / - 1)} . 

In particular, II n := n„ is the space of the unconstrained polynomials of degree at most n. 
Obviously, dim LT^'^ = n — k — I + 1, and the Bernstein polynomials , B% +1 , . . . , -B™_;} 
form a basis of this space. There is a unique dual constrained Bernstein basis of degree n, 

D^ k ' l \x; a, /3), D^ k /\x; a, (3),..., D^ l \x; a, (3), 

satisfying 

Dp Kl \ B%) = Sij = k, k + 1, . . . ,n - I). 

(n k I) 

Recently, we have shown that the polynomials D\ ' can be expressed in terms of dual 
Bernstein polynomials without constraints. Namely, we have the following result. 

Theorem 3.1 ([16J). For i = k, k + 1, . . . ,n — I, the following formula holds: 

(3.1) Z)f ' M) (x;a,/3) = U t ■ x k (l - x) 1 D^~ l (x; a + 21,(3 + 2k), 
where 

I n — k — l\ ( 1 

(3.2) U, : ( . _ k 

Now, using expansion (|2.8p in (|3.ip . the following result is easily obtained. 

Theorem 3.2. The constrained dual basis polynomials have the Bezier- Bernstein represen- 
tation 

n—l 

(3.3) D^ Al) (x;a,/3) = Y,Cv(n,k,l,a,f3)B^(x) (k<i<n-l), 

j=k 

where 

(3.4) C i:j {n,k,l,a,l3) := Ui UjCi_ ktj _ k (n - k - I, a + 21, /3 + 2k), 
notation used being that of (|2.9p and fl3.2[> . 
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... 

Ckk Cfc,fc+1 • • • Cfc,n-i 

Ck+l t k Cfc+l,fc+l • • • Cfc+l 1?1 -Z 



C n -i t k C n _; ; fc + i . . . C n _; jn _; 

... 

Observe that the quantities Cy can be put in a square array (see Table [I]) . 
Combining (|3.4p with Theorem 12. 2 \ we obtain the following theorem. 

Theorem 3.3. Quantities Cij = Cij(n,k,l,a, f3) defined in (|3.4p satisfy the recurrence rela- 
tion 

Ci+ij =^y{(* " j)( 2i + 2j-2n-a + f3) Cy 

(3.5) + 5*(i) Gij-x + A*(i)C iji+1 - B*(i)Ci-i d }, 

where k < i < n — I — 1, k < j < n — I and 

A*(u) := (u - n){u -k + l)(u + k + (3 + l)/(u + 1), 

:= u(u — n — I — a — l)(u — n + I — l)/(u — n — 1). 

PFe adopi t/ie convention that Cij := if i < k, or i > n — I, or j < k, or j > n — I. The 
starting values are 

(3.6) C ki =W-- U l (j = k,k+l,...,n-l), 



where we use notation (13.2D . and 

/ n\-\-l) k (a + 2k + 2l + l) n _*_,(fc + /3 + 2) n _ 



(3.7) W := , , 

K J \kj B(a + 2l + l, l3 + 2k + l){n- k-l)\ 

Using equation (I3.5p . one may obtain the element Cj+ij in terms of four elements from the 
rows number i and % — 1 (see Table [2]). Now, the C-table can be completed very easily in the 
following way. 

Table 2. The cross rule 
Cij-i Cij Cjj_|_i 



BEZIER REPRESENTATION OF THE CONSTRAINED DUAL BERNSTEIN POLYNOMIALS 7 

(n k I) 

Algorithm 3.4 (Computing the Bezier coefficients of the polynomials D- ' ' ). 

(1) Compute quantities Ckk, Cfc,fc+i> . . . , Ck, n ~u filling the first row of the C-table, by the 
formulas 

^nVV n V 1 (-l) n ~ k ~ l (o- + 2k + 2l + !)„_*_, 



(3.8) C, 



k,n—l 



kj \l J B(a + 2l + l, f3 + 2k + l)(n-k-l)V 



r _ U-n)(j-k + l)(j + P + k + 2) 
(3.9) (j = n-l-l,n-l-2, ...,k). 



(2) For i = k,k + 1, . . . ,n — I — 1 and j = k, k + 1, . . . , n — I, compute Cj+ij, using the 
recurrence (13.51). 



Remark 3.5. Obviously, in case k = I = 0, the above scheme allows to compute the Bezier 
coefficients Cij(n, a, f3) of the unconstrained dual polynomial Df(x; a, ft) (cf (|2.8[) ). 

Remark 3.6. If a = (3 and k = I, we have the following symmetry property of the C-table: 

Cij = C n —i jn -j (k < i,j < n — k). 

This can be verified easily, using the definition (|3.4p of Cij as well as the identity 

Q m (n — i; a, a, n) = ( — l) m Q m (i; a, a, n). 

Thus, the cost of computing the C-table can be halved in this case. 

Remark 3.7. Jiittler [6] studied the case of k = I and a = j3 = and gave a rather complex 
recurrence relation in i and k, for the coefficients cfj in 

n—k 

Dp k ' k \x; 0, 0) = B j( x ) (k<i<n-k). 

j=k 

Rababah and Al-Natour [12] extended his approach to the case of arbtitrary a, /3 > — 1. 

4. Applications 

Consider the following approximation problem. Given the function / defined on [0, 1], we 
look for a polynomial P m G LT m , 

m 

(4.1) P m (i):=J> 

which gives minimum value of the squared norm 

(4.2) || / _ jPm ||2 2 = (/ _p m)/ _ j p m) 

with the constraints 

f {i) (0) = P^(0) (i = 0,l,...,k-l), 

(4.3) 

/W(l) = Pi j) (l) (j = 0,l,...,Z-l), 

where k + I < m. 

This general model contains, for instance, 
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(1) multi-degree reduction problem, where / = L n is the polynomial of degree n, n > to, 
given in Bezier form, 



n 

n 
i 

i=Q 



(see, e.g., [16 1 and references given therein); 

(2) computing roots of polynomials by the method of clipping [TOJ O [TO], which uses a 
specific variant of the multi-degree reduction with low value of the degree m, say, less 
than five; 

(3) polynomial approximation of the rational Bezier form, where / = Rn, 



(see, e.g., [9] and references given therein). 

Theorem 4.1. The coefficients po,Pi, ... ,p m of the polynomial (|4.ip minimising the error 
()4.2p with constraints (|4.3|) are given by 

(4.4) pi =^LZ^/(0(o) - £(-1)*+' (%j (< = 0,1 

(4.5) p m _, =(_!)* i!?LZi)l / (i)( 1 ) _ J2(-iy ( l ) Pm - i+j (i = 0, 1, . . . , l - 1); 

m— 2 / fe— 1 m \ 

j=fc \i=0 j=m-l+lj 

(4.6) (i = jfe,Jfe + l,... 1 m-Z), 
where Cij(m,k,l,a, f3) are introduced in (]3.3[) anc? 

'm\fmX~ (-iy~ k (k - j) m _ fc _ i+ i 



(4.7) 



1 jV V * -/ (i — — A;)!(m — Z — i)! 
(a + / + l) m _ j (/3 + fc + l)j 



(a + i + lV-iC/S + fc + lJi" 
Proof. Recall that (see, e.g., [U p. 49]) 

i 



^(o)^E(-^(i)-. 



Using the above equations in (|4.3p . we obtain the forms (|4.4|) and (|4.5h for the coefficients 
po, pi, • • • ,Pfc-i and p m _/+i, . . . ,p m -i,Pm, respectively. 
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The remaining coefficients p, are to be determined so that 

ra—l 

11/ -Pm\\l 2 = \\W-Y,PiBT\\l 2 

i=k 

has the least value, where 

(k-l m \ 

E+ E U*r- 
j=0 j=m-l+lj 

Remembering that B™ and £)( m ' k ' 1 ^ (k < i < m — I) are dual bases in the space 11^'' , we 
obtain the formula 

m—l I k— 1 m \ 

Pi = (W, D^' k ' l) ) = Cij(m, k, I, a, /?) {/, Bf ) - £ + £ p^-, 
j=fc \i=Q j=m-l+lj 



where 



ify := (Bf , .Of*"' 



Formula (14. 7|) was obtained in [9]. Equation (14.60 now readily follows. □ 

Remark 4.2. Evaluation of quantities (/, -B™), which are present in equation (|4.6p . is a simple 
task in case of multi-degree reduction or computing zeros by Bezier clipping, as we have 



(/' B T) = < L «' B T) = E ^ W*. 

and 



i 

i=0 



<BT.^)-^ + i./« + i)0(7) ( ° + ^^ 

In case o/ approximation of the rational function, we have 



(f,Bj l ) = (R n ,B^) = Y / ^r i I, 



i=0 



where the integrals Jjj := -B™ / Yl^o^h B h) s ^ ou ^d be computed using numerical meth- 

ods (see, e.g., [9],). 



Conclusions 

We gave explicit formulae for the Bezier coefficients of the constrained dual Bernstein 
basis polynomials in terms of the Hahn orthogonal polynomials. Using difference equation 
satisfied by the latter polynomials, we obtained a recursive algorithm to compute efficiently 
these coefficients. Applications of this result to several problems is discussed, namely to the 
multi-degree reduction problem for the Bezier curves, computing roots of polynomials by the 
method of Bezier clipping and approximation of the rational Bezier curves by Bezier curves. 
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